Non-Ohmic hopping transport in a-YSi: from isotropic to directed percolation 
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Electrical transport has been investigated in amorphous Yo.19Sio.8i, from 30 mK to room temperature. 
Below 2 K, the conductance G exhibits Efros-Shklovskii behavior G ~ exp[— (Tq/T) 1 / 2 ] at zero electric field, 
where conduction is expected to occur along a very sinuous path (isotropic percolation). The non-linear 
/-^characteristics are systematically studied up to very high fields, for which the conductance no longer 
depends on T and for which the current paths are expected to be almost straight (directed percolation) . We 
show that the contributions of electronic and sample heating to those non-linearities are negligible. Then, 
we show that the conductance dependence as a function of low electric fields (F/T < 5000 V m _1 K _1 ) is 
given by G(F, T) = G(0,T) x exp(-eFL/fc B T). The order of magnitude (5-10 nm) and the T dependence 
(~ T -1 / 2 ) of L agrees with theoretical predictions. From the To value and the length characterizing the 
intermediate field regime, we extract an estimate of the dielectric constant of our system. The very high 
electric field data do not agree with the prediction 1(F) ~ exp[— (Fq/F) 1 '] with 7/ = 1/2 : we find a F 
dependence of 7/ which could be partly due to tunneling across the mobility edge. In the intermediate 
electric field domain, we claim that our data evidence both the enhancement of the hopping probability 
with the field and the influence of the straightening of the paths. The latter effect is due to the gradual 
transition from isotropic to directed percolation and depends essentially on the statistical properties of the 
"returns" i.e. of the segments of the paths where the current flows against the electrical force. The critical 
exponent of this returns contribution which, up to now was unknown both theoretically and experimentally, 
is found to be (3 = 1.15 ± 0.10. An estimation of the returns length is also given. 
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I. INTRODUCTION 



The electrical transport in disordered insulators such as doped semiconductors or amorphous systems, 
has been widely investigated for several decades as a convenient way to probe the properties of Anderson 
insulators [Q . At finite temperature T, the conduction results a priori from all the electronic hops between 
any localized states i and j separated by a distance and a difference in energy £Jy. Starting from 
the probability of such a hop pij oc exp — (ry /£ + Eij /ksT) (with £ the localization length and ks the 
Boltzmann constant), Mott Q introduced the key concept of variable-range hopping (VRH) stating that 
transport is dominated by the most probable hops whose characteristic length r m oc T _1 ^ d+1 ^ and energy 
E m oc T d /( d+1 "> depend on T. As a result, the conductance G is predicted to follow 

G = G cxp[-(^) 7 ] (1) 

where Go and To are material dependent constants, and 7 = l/(d+ 1) depends on the dimensionality d 
of the sample. Later, the percolation theory has been used to derive more rigorously Eq. (1) || Q. In 
addition Efros and Shklovskii || found that electron interactions result in a "soft" Coulomb gap at the 
Fermi level, leading to 7 = 1/2. It is currently admitted now ||] that Mott and Efros- Shklovskii laws 
are limiting cases of VRH. Both of them were observed experimentally as their relevance depends on the 
relative magnitudes of the hopping energy and Coulomb gap. 

A. Electric field effects in hopping transport 

The case of very high fields F, eF^/ (kBT) 3> 1 (e = electron charge) was studied through differ- 
ent theoretical approaches @]-JlJ] which all predict that the current / should no longer depend on the 
temperature (activationless hopping) and behaves according to 



I = Io exp 



(2) 



FJ . 

where Io and F are constants and 7/ = 7 (see Eq. (1)). In spite of their similarity, Eqs. (1) and (2) 
correspond to completely different current path topologies. Eq. (2) is obtained by considering the shortest 
hops along which the electrostatic energy gain overcomes the energy fluctuations due to disorder |l3|, |7|. 
Thus in the very high field case, the hops occur in the direction of the field (directed percolation). On the 
contrary, for F — ► 0, the current flows through a network of random impedances along very sinuous paths 
whose returns and meanders are such that any hop more resistive than the VRH prediction is forbidden. 
A quantitative description of these paths can be obtained by using the isotropic percolation theory fbl 
II- 

The intermediate field case eF£/ (ksT) < 1, where neither T nor F can be neglected is the most difficult 
one to handle theoretically, and the various models lead to different predictions that we summarize in the 
general equation giving the current 

«™-'-- Hir M^SMari} (3, 

where A, B, a, af and (3 are parameters, A and B being positive. The first term A(eFt;) a / '(feT)" 
expresses the enhancement of the hopping probability when F grows. Depending on the authors, a and 
ai range from a = a' = l/(v + 1) Jl0| | (where v ~ 0.88 for d = 3 is the critical index of the correlation 
radius) to a = a! = 1 ||], p3| , and even a = 2 and a' = 9/4 for d — 3 ||, [Q. The discrepancies between 
these results come from the fact that the net current hj between two sites may increase exponentially 
with F or be insensitive to F, depending on the relative positions of the site energies. The second term 
B(eF^/{kBT))^ was predicted only by Bottger, Bryksin et al. who addressed the additional 

problem of the "returns" , i.e. of the segments of the paths where the current flows against the electrical 
force. Their model states that the length of these returns decreases gradually when F increases. We 
show in this paper that our data evidence the existence of both terms of Eq. (3) with a = a' = 1 and 
P = 1.15 ± 0.10. 
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B. Theories of hopping transport at intermediate fields 



We recall here some theoretical points which are important for the interpretation of our data. The 
very sinuous current paths at F — > can be pictured by using the "nodes- links-blobs" concept Q |l5| , 
derived from the isotropic percolation theory : at small scales the current path consists of a fractal-like 
network of connected conductances up to the "blob" scale C p , and the blobs are the links of a homogeneous 
network. In the linear regime, the conductances inside a blob are much larger than the key conductances 
Gkey ~ exp — (T /T) 7 lying at the end of each blob. Gk ey sets the overall sample conductivity since at 
scales larger than C p the system is homogeneous. For the intermediate field case, where neither T nor F 
can be neglected, three main effects have to be considered : 

(i) for a site i belonging to the percolation path, the density of sites j yielding a non negligible con- 
ductance dj is modified by F; 

(ii) a field-induced rearrangement of the charge along the current paths can take place : the local 
chemical potential /Ltj can fluctuate widely around the mean value given by eFrf, 

(iii) the lengths of the "returns" , i.e. the segments where the current has to go against the electrical 
force, decrease when F increases. 

Most authors disregarded (iii) and focused upon (i) and (ii), yielding predictions summarized by Eq. 
(3) with B = 0. However, they gave diverging predictions depending on the relative influence given to (i) 
and (ii). For example, Pollak and Riess || find a = a' = 1 and A = 3/32 in a model where (i) plays the 
central role, and the correlations between nearest neighbors conductances are taken into account. They 
find that (ii) is quite negligible (Sfii < UbT). This contradicts Shklovskii model |L0| [A ~ 1 and a = a' = 
1/(1 + 1/) ~ 0.53 for d = 3) which emphasizes (ii) (SfXi 3> fesT) by concentrating almost all the electrical 
potential drops existing at the blob scale C p upon the key conductance. The disagreement between these 
predictions was never completely explained |is|| even if numerical simulations by Levin and Shklovskii |20|] 
suggested that the Pollak and Riess result holds in usual experimental conditions (To/T) 1 ' ( d+1 ) < 70 while 
Shklovskii result would hold in the opposite case, i.e. at exceedingly low temperatures. Two other groups 
found different results : Apsley et al. § and Van der Meer Q obtained a = 2, al = 9/4, but their 
methods were criticized since the former did not take into account the directed percolation requirement, 
while the latter uses dimensional "invariants" which were shown later to be F dependent Ell . 

Taking (iii) into account jjl(||2l]]|^] prevents the current from increasing too fast : the corresponding 
second term B^F^/ksT) 13 in Eq. (3) reduces the exponential increase of I due to — A(eF£) a / (fceT) a/ . 
The physical reason is that the returns act as "bottlenecks" for the current. The problem being exceedingly 
complicated, the precise values of a, (3, A, B were found to depend on the various methods used by Bottger, 
Bryksin et al. The last and more accurate results were obtained numerically p^||Q |Q : a = al = 1, 
A = 1/6 , B ~ 0.02, 1//3 ~ 1.1 (this latter parameter is the critical exponent of the returns length A, see 
below Eq. (12)). The theoretical value of (3 is given within an error of 50% due to numerical uncertainties. 
Yet, the precise value of /3 is of major importance since for (3 < a, the differential conductivity a — dl/dV 
decreases at low fields to reach an absolute minimum at a finite F. On the contrary, in the case (3 > a, a 
increases at low fields and has no absolute minimum at finite field. 

The above mentioned theories disregard the possible independence of the carrier density due to trapping 
in dead ends of paths, as well as the possible -F-dependence of the localization length £. Fortunately, both 
effects should be negligible in our case. Indeed, statistically as many holes as electrons become trapped 
when F increases, as in our case the density of states is symmetric around the Fermi level (parabolic 
Coulomb gap) (23|. The £(F) dependence should be negligible Jgjj] as long as Ep 3> eF£ where Ep is the 



kinetic energy at the Fermi level: this is the case in our sample [25 



C. The experimental situation 

On the experimental side, I-V non linearities were investigated both in amorphous materials |2(|-[|34| 
and in doped crystalline semiconductors p5j-p7|. In most of these works, the authors focussed either on 
very high fields or on intermediate fields. For intermediate fields, the data Q Q J^H j}5)-|42 IHl were 
analyzed using the electric field dependence predicted by Hill jl3| , Pollak and Riess M, and Shklovskii fl(J 

In (G(F, r)) = In (G(0, T)) + f± = In (G ) - @) " + f± (4) 
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where L is a length related to the hopping distance r m = (£/2)(To/T) 7 . Hill[T^] and Pollak and Riess|| 
predict L = C.(£/2)(T /T)t with C = 0.75 or 0.18, while Shklovskii predicts L oc T^ 1 ^ |§ with 
v = 0.88. Eq. (4) results from Eq. (3) by neglecting the returns term (B = 0) and by assuming that 
the F and T dependences of G and / are close to each other. The experimental values of L and their T 
dependence often disagree with the theoretical predictions [^6|-|3^|pt|. The claim by Wang et al. Q 
that the non-linearities in Ge:Ga were due to carriers heating ]64| rather than to the preferential hopping 
along the electric field which leads to Eqs. (3) or (4) triggered hot-electron analyses of the data p7|-|53|. 
They show that the relative contribution of carriers heating and electric field effects has to be considered 
at very low temperature. In very high electric fields, the expected activationless conduction (see Eq. (2)) 



has been observed by several authors [|30|-|32 |12| [j43| J46|]. The equality of the VRH and activationless 
exponents (7 and 7/ in Eqs. (1) and (2)) has been checked for 7 = 1/2 in several systems j3^] [ p3[ , but 
7/ = 1/4 and 7 = 1/2 has been found in amorphous GeCupl[|. 

In this work we systematically investigate the whole set of non linearities ranging from the linear VRH 
regime at F — > 0, to the activationless conduction at very high fields : in our low temperature domain 
(0.03 K < T < 1.3 K), F/T varies from 1.5 x 10 2 Vm-'K" 1 to 1.4 x 10 6 Vm-'F 1 . Our samples 
are made of amorphous Yo.19Sio.8i- In comparison to most of the experimental works performed in this 
low temperature range ( < 1 K) |37j p3[-|47], wc go to higher values of F and F/T, thus allowing an 
investigation of both the intermediate and the very high field regimes. At very low temperatures, this has 
been done by Rosenbaum et al. on Si:P p7j| , however the authors found an abnormally large value of L. 
At higher temperatures, these two regimes were investigated for a-GeCu at 1.6 K < T < 110 K |n] and 
ZnSe at 1.6 K < T < 4.2 K p2[ . The first paper raised the question of unexpected exponents values (in 
the L vs. T and ln(G) vs. F dependences) while the second one found a puzzling decrease of L when F 
increases. Clearly, new data are needed. In this work, we test the validity of Eqs. (l)-(4) and show that 
they are in good agreement with our data for intermediate fields. In particular we find L ~ C.r m @§. 
We also check that the contribution of heating effects to the non-linearities is negligible. Our main result 
is that Eq. (3) with a = al is more suitable than Eq. (4) to fit the data : this allows us to estimate the 
size of the returns and the value of the related critical exponent /3[|l6|-[|l8|. Finally, we show that in the 
activationless region, the I vs F dependence is more intricate than the law given by Eq. (2). 

The paper is organized as follows. In section II we describe the experimental setup and method. The 
conductance as a function of the temperature in the region F — > is presented in section III. In section IV, 
we consider the very high field case. Having extracted the physics of the two extreme F regimes, we finally 
turn in section V to the intermediate field case. Finally, we summarize our main results in section VI. In 
appendix A we show that heating mechanisms cannot explain the I — V non-linearities, while appendix B 
is devoted to the extraction of the localization length from the data. 

II. EXPERIMENT 



As shown in Fig. 1, each amorphous Y K Sii_ K sample studied (x ~ 0.19) is deposited on a sapphire 
substrate on which two interdigited gold electrodes were evaporated and etched. The spacing between the 
electrodes is I = 128 /mi. The Y x Sii_ x layers were obtained by Argon plasma sputtering on the sapphire 
substrate (cooled at 77 K to prevent Y or Si aggregates formation), using a YzSii-z source. The thickness 
of the layers is 9 /tm, much larger than any hopping length, hence electrical transport is three dimensional. 
Previous studies of electrical transport in such Y^Sii-^ samples |54|-|39|], have shown that a small variation 
of x allows to cross the metal-insulator transition : above x ~ 0.22 the samples are metallic; while for lower 
values of x, the samples exhibit an insulating behavior, i.e. a divergence of the resistance for T — + 0, in 
agreement with the VRH predictions j55|-|58|, except for the very weakly insulating samples at the lowest 
temperatures |59| . In the present work x ~ 0.19 and the samples are strongly insulating, exhibiting VRH 
below T ~ 2 K. Since the two samples wc used behaved similarly, we report only the data of one of them. 

As shown in Fig. 1, a great attention was paid to the thermalization of the sample because of a possible 
heating interpretation of the I-V curves non-linearities. As explained in Appendix A, we measured the 
thermal conductances involved in our sample and concluded that heating effects were irrelevant in our 
sample: neither heating of the whole sample nor carrier heating can explain the I- Vnon linearities reported 
here. 
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Electrical contacts were made by ultrasonic soldering of two gold wires on the evaporated electrodes. 
The data were obtained by setting a given voltage V — Fl on the sample, and by waiting for a time r 
before measuring the current /. We paid the greatest attention to choosing r large enough to let the 
current settle, e.g. we took r(T < 800 mK) > 100 s. The temperature was registered for each /(V) point, 
and its stability proved to be better than 1%. The coaxial cables used were previously tested alone (in 
an open circuit, at all temperatures) to check that their leakage current was negligible. Let us note that 
it is unlikely that the gold-YSi electrical contact resistance plays a role since the contact area is very 
large and in previous similar works |58[ |, it was shown to play no role. The fact that the temperature 
dependence G(0, T) we observe follows a VRH law (see next section) is an additional indication that the 
contact resistance is negligible. Finally, the symmetry of the I- V curves with respect to current and voltage 
reversing was checked. 

III. BEHAVIOR ATF^O 

The /-^characteristics are shown in Fig. 2 where, for clarity, we report only the curves for 18 different 
temperatures out of the 36 ones we measured. As shown in the inset, the resistance (R) - temperature law 
is simply activated from room temperature to ~ 10 K : we find R = iiiexp(— T%/T) with T\ — 40 ± 4K 
and R\ ~ 88 ± 2 f2. For T < 8 K the divergence of the resistance when T — > becomes weaker (see inset 
of Fig. 2), which is a standard indication that VRH takes place at low temperature ||. Let us note that 
simple activation for T > 10 K is usually interpreted as due to thermal excitation from the Fermi level Ep 
to the mobility edge Eyiob where electronic states are delocalized; thus fcgTi = E^ob ~~ Ep. 

To study the temperature dependence of the conductance at F — > 0, we plotted our data in the 
{F/T, \n(V/I)) plane (see Fig. 3). The G(F — > 0) values were obtained from an extrapolation of the 
curves towards low fields. To allow for a precise extrapolation, this procedure was restricted to the 
0.4 K < T < 1.3 K cases. The extrapolated values were obtained from a linear fit of the first points in 
each curve (lower F/T values). We checked that those values did not depend significantly on the number 
of points selected in the fit. The expected VRH temperature-conductance law is given by Eq. (1) where 
7 = 1/2 or 7 = 1/4. Discriminating between these exponents is hardly achieved by looking for a straight 
line in the bx(G(F — > 0)) versus T ' plot, as the lines always seem straight. The comparison between the 
different exponents 7 is improved by plotting a normalized value of T~ 7 , i.e. (T~ 7 — T^jJ n )/ {T^p ax — T^jJ n ) 
where Tuin and Tuax are the two extreme values of T selected for the plot. In Fig. 4, we see that 7=1/2 
is favored in comparison with 7 = 1/4. A fit of \n(R(F — > 0)) vs. T using Eq. (1) in which Go and To 
are free parameters yields a residue (measured minus fitted value) shown in the inset of Fig. 4 : clearly, 
there is a correlation between the residue and T only for 7 = 1/4. The normalized \ 2 P er point is 1.1 for 
7 = 1/2, and 9.0 for 7 = 1/4. If 7 is a parameter of the fit, we find 7 = 0.47 ± .02. We thus conclude that 
7 = 1/2, and using this value in the fit, we get T = 257 ± 1.5 K and R = 1/Gq = 92 ± 5 £1. Following 
Efros and Shklovskii, 

2 9e 2 

k B T = — (5) 

where eoe r is the dielectric constant of the system. We shall see in Appendix B that the non-linearities 
analysis together with Eq. (5) allow to extract an information on £ and e r separately. 

IV. RESULTS AT VERY HIGH FIELDS 

In the very high fields case, eF£/ (fcgT) S> 1, the theory states that the current results from hops 
between sites i and j, such that E^j = eFrij. For these activationless hops, oc exp (— 2r , i 7 /£), which, 
using a maximization procedure leads to Eq. (2) with0 

F = 6 

where N% is a numerical constant which can be calculated within the directed percolation theory. For 
Mott's VRH (7 = 1/4) Pollak and Riess obtain N 2 = 4.8. For the Efros-Shlovskii VRH (7 = 1/2), N 2 has 
not been calculated to our knowledge. 
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Our lowest temperature data should belong to the very high field domain, as the two I- V curves 
corresponding to T = 29 mK and T = 49 mK are identical. Indeed, we shall see that the low and 
intermediate field data analysis lead to £ w 7 nm (see Appendix B), hence eF£/(ksT) > 50 for T — 29 mK. 
The 29 mK data can be rather well fitted by using Eq. (2) with 7/ = 1/2, Iq and Fq being free parameters. 
It yields F = (3.8 ± .05) x 10 7 V/m, hence 4 < N 2 < 11 by using the upper and lower values of £. The 
normalized \ 2 P er point is however 290. Thus, in spite of the reasonable extracted Fo and iVjj values, the 
relevance of the fit must be questioned. Indeed, if we fit the whole I(F, 29 mK) data using Eq. (2), 7/ 
being a free parameter, we get 7/ = 0.65 ± 0.01 and a x 2 per point of 94. 

To investigate the problem, we give in Fig. 5 the F dependence of the effective exponent 7 e // defined 

by 

qi ( dlnl \ 

7e// " dbxF ■ [7) 
If 1(F) is given by Eq. (2) , then j e f f = 7/ . For T > 200 m,K, j e f f is always negative, while for T < 200 mK, 
j e ff it is negative for the lower F values. This can be related to the fact that in the intermediate field 
region, the 1(F) dependence should obey Eq. (3) for which the expected j e ff is F-dependent. For the 
lowest temperatures, 7e// becomes positive when F increases, but there is no saturation at 7e// = 0.5. 
Instead, a maximum is reached at a value between 0.5 and 1, even if a convergence of Jeff towards 1/2 
might exist at fields larger than 4 x 10 4 V/m. 

We suggest that these values of 7 e // between 0.5 and 1 are due to the fact that the very high field 
conduction results from an interplay between activationless hopping (for which J e ff = 1/2) and tunneling 
across the mobility edge. Indeed, the latter process leads to the simple activated law |3q| , |30|j 



exp 



J_ /2m( fc .TQ 
V 3eF 7 V h 2 y 







= exp 


. F . 



where fcsTi is the energy difference between the Fermi level and the mobility edge. Since T\ ~ 40 K (see 
Sec. Ill), F x = 1.4 x 10 6 V/m. The fit of the 7(F,29mK) data for the interval F > 2.8 x 10 4 V/m where 
j e ff > 0.5, using I — I\ exp (—F\/F), yields a normalized x 2 P er point of 40 and F\ = (5.5±0.05) 10 5 V/m, 
which is only a factor 2.5 below the calculated value. We cannot interpret the j e f / (F) curves more precisely 
since, to our knowledge, there exists no theory taking into account both the disorder in the localized band 
states and the tunneling through the Mobility Edge. 

It is interesting to note that if activationless hopping was the only transport mechanism in the very 
high field regime, the critical value (F/T) c at which the transition from Eq. (3) to Eq. (2) occurs 
would be temperature independent : this can be readily seen by equating the rhs of Eqs. (2) and (3) at 
F/T — (F/T) c and disregarding pre-exponential factors. If we define, for any given T the experimental 
critical field value F c by I(F C ,T) = kI(F c , 29 mK), with k = 1.1, we find that in our sample, (F/T) c is 
clearly T-dependcnt: it increases by a factor 3.5 when the temperature decreases from 413 mK to 74 mK. 

V. RESULTS AT INTERMEDIATE FIELDS 

We can see on Fig. 3 that our experimental points line up in almost parallel curves in the (In (i?), F/T) 
plane. Clearly, they are not really straight lines, hence Eq. (4) does not hold with precision. However, in 
the spirit of the majority of previous works, we first use Eq.(4) and extract the length L(T) defined by 
In [G(F)/G(F — > 0)] = eFL(T) / '(fcgT) for F just above the linear regime. We show in a second step that 
Eq. (3) is much more relevant to account for our data. 



A. Extraction of L(T) using Eq. (4) 

The slopes of the In (R) vs. F/T curves for F — *■ were extracted for the different temperatures where 
it was possible, i.e. 0.4 K < T < 1.3 K (see Fig. 3). A linear fit yielded a normalized % 2 per point close 
to 1 only for F/T < 5000 Vm _1 K _1 . We restricted the L extraction to this interval where the curves are 
linear. The result is shown in Fig. 6. A fit of these points using the law L = Lq (Tq/T)^ where Lq and £ 
are free parameters gives C — 0.65 with a normalized \ 2 P er point of 0.7. Constraining the fit with ( = 0.5 
leads to Lq = (0.52 ± 0.01) nm and a normalized x 2 P er point of 1.1 : the resulting curve is shown in 
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Fig. 6. As a consequence, we can state that the Shklovskii prediction |l0| ( — 2-f — 1 does not hold for 
our system, while the predictions of Hill [jl3| and Pollak and Riess ||, L oc r m leading to C = 7 = 1/2, are 
in agreement with our results. It is interesting to note that the values of L(T) remain roughly the same if 
they are extracted from linear fits of the whole ln(i?) vs. F/T curves, although these fits yield a normalized 
X 2 per point varying between 3 and 250. The localization length £ extracted from L = C.(£/2) with the 
Pollak and Riess value C = 0.17 || yields £ ~ 6nm. We shall come to this point more in detail later. 

B. Extraction of the parameters of Eq. (3) 

We turn to the main result of our study, namely the relevance of the Bottger et al. prediction |l7| ] 
concerning the "returns" and summarized in Eq. (3). 

B.l. F/T is the unique relevant parameter 
We define, so as to test Eq. (3), 

*(F,T) = Mki( I{F ;l\ ] ) ,with J £ . s .(F,T) = ^exp[-(^) 1 /2] (9 ) 
V Jo \ Ie.s. \-F , J ) J H 1 

where V = Fl, R a = 92 ft and T = 257 K (see Sec. HI). It is easy to show that if Eq. (3) holds in what 
concerns the G(T, F) dependence (instead of I(T, F)), * is just equal to A(eF£) a /{k B T) al '~B(eF^/k B Tf '. 
We have chosen to perform our analysis on G(T,F) = I(T, F)/ (Fl) rather than on I(T,F) because all 
the above depicted theories do not allow to recover a linear regime when F — > : in this case, Eq. (3) 
gives I(F — > 0) = I\ exp [— (T /T) 7 ] which is not consistent with Eq. (1). Note that this difficulty has 
been ignored in almost all experimental works dealing with the intermediate fields case, as they focus 
upon G(F) rather than on 1(F). Numerically, our choice is justified by the fact that the exponential 1(F) 
dependence in Eq. (3) is much faster than the contribution of the denominator of G(F) ~ 1(F) /F (see 
also |6l|). 

To perform our analysis of the intermediate field case, we suppressed the data points corresponding to 
the very high field case (see Sec. IV). This was done by selecting 7164 I(F, T) points among the whole set 
of 8189 points, on the criterion I(F, T)/Ivhf(F) — 1 > 10%, where Ivhf(F) is the very high field current 
defined by Ivhf(F) = I(F,T = 29mK). The 10% criterion is rather arbitrary, but we checked that our 
results do not depend on its precise value, provided it remains above 5%. Then, \I> is calculated for each 
of these 7164 points. As shown on Fig. 7, ^(F, T) is a universal function of F/T over more than 3 orders 
of magnitude, while the 7164 experimental points correspond to 36 different temperatures ranging from 
49 mK to 1.29 K. 

The fact that ^(F, T) depends only on F/T implies a = a' , which definitely shows that both Apsley § 
and Van der Meer |12| models are not relevant for our sample. Moreover, the log-log slope of ^(F/T) at the 
lowest F/T values is very close to 1 and clearly larger than \/(v + 1) ~ 0.53 which means that Shklovskii 
prediction |l(| does not hold. This is not surprising since in our experiment (Xb/T) 1 / 2 < 70 which is a 
domain where the numerical work of Levin and Shklovskii |20| leads to a. = 1. However our data are poorly 
fitted by the Hill Jl3| and the Pollak and Riess j§ prediction * - F/T since the log- log slope of ^(F/T) 
decreases as F/T increases. More generally, a fit of ^(F/T) with a unique power law (F/T) a gives a poor 
agreement whatever a", especially because of the points in the region F/T > 3 x 10 4 V K^m- 1 . The only 
remaining prediction is thus the one of Bottger et al. where a — 1 and (3 is unknown theoretically. In the 
next section this prediction will be assumed and the best value of (3 will be sought. Quite interestingly, 
the fit of ^ we obtain in the next section is, by far, much better than all the "reasonably simple" fits we 
tried: for example attempting a polynomial fit with ^(F/T) 1 gives a poor agreement even if m is as 
large as 5. 

B.2. Extraction of the critical exponent [3 assuming a = 1. 
We thus have to fit the ^(F/T) points using the function 

r£-A(£) fl (io) 
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with r, A, as free parameters. Since the ^(F/T) values extends over several orders of magnitude, we 
perform the fit on In ($), i.e. we request that a relative error £ per point is minimal, with £ defined by 



£ 2 = 



+x 



hi 



7164s 2 
r(F/T) J -*((F/T) 



J= i,7i64 ^(F/T), - A(F/T)^ 



A(F/T)f 



(11) 



where x is a weight to be chosen in the (0, 1) interval and s 2 — ((ln(5')— < ln(^) >) 2 ) — (0.08) 2 is the 
variance of ln(\I/) drawn from the data. Note that the first term in the rhs of Eq. (11) is the relative error 
per point on ^ while the second one is the relative error per point on the "Bottger second term" of Eq. 
(3), A(F/T)@ . The minimization of those two terms has opposite effects on the optimal value of (3. Indeed 
keeping only the first one (i.e. setting x = 0) leads to (3 = 1.38 ±0.03 and yields a good agreement between 
the ty((F/T)j) points and the fitted curve, but leads to a discrepancy between T(F/T)j — 5 , ((F/r) J ) and 
A(F/T)j which becomes clearly too large when F/T < 5 x 10 4 V K _1 m -1 . Conversely, setting x = 1 
leads to f3 = 1.03 ± 0.03, with T(F/T)j - V((F/T)j) ~ A(F/T)f, but a discrepancy appears between the 
^((F/T)^ points and the fitted curve, especially at high F/T values. Keeping x in the (0.05, 0.95) interval 
leads to optimal values (3(x) in the (1.05, 1.15) range, i.e. to an uncertainty on (3 not very much larger 
than the one obtained for a given x. We finally keep the larger error bar on (3 and find : (3 = 1.15 ± 0.1, 
T = (1.25 ±0.5) x 10~ 5 K m/V and A = 1.4 x 10- 6 (m K/V) 13 (A^/r = 0.70 ±0.07). The solid line on 
Figs. 7 and 8 is the fit resulting from (j3 = 1.15; V = 1.25 x 10~ 5 KmV _1 ; A = 1.4 x 10" 6 (Km V~ 1 ) /3 ) : 
the agreement with the data is good (£{x — 0.5) = 0.72) and the only systematic deviation occurs at the 
very few highest F/T values where the transition to the very high fields regime occurs. 

The comparison of Eqs. (3) and (10) gives T = Ae^/ks- From Y = 1.25 x 10 -5 m K/V we deduce 
= 1.1 nm. Unfortunately, the precise value of A is not known in our case of Efros Shklovskii's hopping 
(7 = 1/2). As shown in Appendix B, using To = 257K, A^ = 1.1 nm as well as considerations about the 
dielectric constant c r we conclude that £ as 7 11111. We cannot be more precise due to the various unknown 
numerical factors involved in the predictions we used, but we note that these values of £ compare favorably 
to previous results obtained on Y x Si\- x samples much closer to the metal-insulator transition |fi8|| , |59|| 
where larger values of £ were found (a few tens of nm) . 

We compare now the analysis performed just above with the one carried out using Eq. (4) which 
yielded L = L (To/T) 1 / 2 with L = 0.52 nm, and ?~6nm, using L = C (£/2) and C = 0.18 (see Sec. 
B.). Neglecting the weak difference between the G(T, F) and the I(T, F) dependences, the comparison of 
Eqs. (3) and (4) leads to Lq = A£ if the additional term B{eF^J '(ksT)) 13 in Eq. (3) is assumed not to 
change Lq. Clearly, we have a discrepancy by a factor 2 between Lq and At; = 1.1 nm : this was checked 
to come from the B{eF^/ {ksT))^ term in Eq. (3). Surprisingly both analysis yield roughly the same 
£ « 6 - 7nm for A = 1/6 and C = 0.18. 

B.3. Extraction of the return length A. 

We now turn to the analysis of the "returns" size. As explained above, (3 and A give an information 
about the typical length A of the returns which decreases gradually as F increases. Indeed, according to 
Bottger et al. [^2|, if p = In [I(F, T)/I{\ (p is the argument of the exponential in Eq. (3)), A is given by 

A(P)-2^^W^Y (12) 
\p-p c J V p c ) 

where /1 = 1.0 ± 0.3 and 5 ~ 0.25 are numerical constants |2^], p m is the directed percolation threshold, 
and p c = (To/T) 7 . In the linear regime it is found that p(F — > 0) ~ p c + 1 JlJ], yielding a finite value of 
A(F — > 0). Increasing F up to the intermediate field region leads to a p increase, thus to a A decrease. 
Increasing F further, p reaches p m and A vanishes indicating that the current paths are directed. The 
maximum length A max of the returns is thus obtained just at the onset of the non linear regime, where we 
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can put in Eq. (12) p = (To/T) 7 + 1 and (p m — p)/p c — (p c — Pm)/Pc — N4 with A/4 ~ 0.07 as predicted 
in p2[| . Hence we get 



T e V T 



where the second equality is obtained by using the Bottger et al. prediction A = (N^Se^/ks) 13 §§]. Using 
Eq. (15), we can compare A max to the blob length C p = ^{T /T)^ ( - 1+ ^ : 



A max AH B {ToY^-^ A?k E 



(14) 



where the second equality comes from the fact that in our case, due to the experimental value of (3, A max 
and C p have critical exponents very close to each other since 1/(3 — v is in the (-0.08, 0.07) interval. We 
thus find that A max (T) / C P (T) does not depend on the temperature. By inserting fce/(e£) = A/T in Eq. 
(15), we find A max (T)/C p (T) = AA? /T = (0.70 ± 0.07)A , using our fit result A?/r = 0.70 ± 0.07. 
Hence, A max (T) / C P (T) is not much smaller than 1 which is the highest possible theoretical value. We 
thus conclude that in our sample the importance of the returns is strong. 

This importance of the returns must also play a role on the onset of the non linear behavior which 
occurs at a field F[ im (T). Indeed, according to Bottger et al. Fi im /T = kg/ 'eA max p2| . According to 
Pollak et al., Fu m /T — ks/er m with r m = (£/2)(Xo/T) 7 the hopping length. In our case, this latter 
prediction amounts to Fu m /T > 1000 V Km -1 for T > 500 mK (for lower temperatures we cannot go to 
F — > because of too low currents). We clearly see on Figs. 7 and 8 that the non linear fit obtained above 
extends down to F/T ~ 150 V K _1 m -1 . This is one order of magnitude smaller than the Pollak et al. 
prediction and in quite good agreement with the Bottger et al. result. We thus conclude that the Fu m 
behavior we find is consistent with the above derived results on the importance of returns. 

VI. CONCLUSION 



Our study of the electric field effects in variable-range hopping transport for amorphous Yo.19Sio.8i 
below 2 K exhibits several important features. First, we find that the length L characterizing the interme- 
diate field regime has the order of magnitude and T dependence (~ T -7 , with 7 = 1/2) which is expected 
in the VRH models of Hill |l3| or Pollak and Riess § stating that L ~ r m . Even analyzing our data in the 
framework of the Bottger, Bryksin et al. predictions (Eq. (3)) this result remains true, as the addition of 
the B(eF£/ (ksT))P term in Eq. (3) does not change the order of magnitude of this characteristic length. 
Second, our most important result is that Eq. (3) is much more relevant than Eq. (4) to analyze our 
data, and this shows the importance of the "returns" in the percolation paths of VRH. Furthermore, we 
were able to extract an information on the length of these returns from our experimental results. Indeed, 
our data indicate that the critical exponent of the returns contribution 1//3 is very close to the one of the 
blob length : the returns represent an appreciable T-independent fraction of the whole percolation paths 
lengths at intermediate fields. Third, our very high field data do not follow the expected activationless law 
given by Eq. (2). This could be due to the onset of tunneling across the mobility edge whose interplay 
with activationless hopping has not been theoretically studied yet. 
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APPENDIX A: IRRELEVANCE OF THE SAMPLE OR CARRIERS HEATING 

MECHANISM 



In this Appendix, we first describe the experimental setup allowing thermal conductances measurements 
and then show that heating is irrelevant to account for our I-Vnon linearities. 

The sapphire substrate was hold between four sapphire balls which represent a negligible thermal 
conductance. The thermalization was realised by ultrasonically bonding 43 gold wires (50 (im in diameter, 
a few millimeters long) between the cryostat ("Heat sink" on Fig. 1) and a 7mm 2 gold pad schematized 
on the left of the substrate on Fig. 1. The resulting calculated thermal resistance between the gold 
pad and the cryostat is smaller than 10 4 K/W at T = 50mK and 10 3 K/W at T = 500 mK, leading to 
a temperature difference lower than 0.5 mK between those two points in our experiments. To measure 
the thermal resistance between the cryostat and the sapphire substrate or the sample, we used the 1 M.Q 
heating resistor chip shown on the right up corner of Fig. 1. This 0.01mm 3 device is held only by two 
Al wires (25 /im in diameter, 5 mm long) which are superconducting below ~ 1 K. At low temperature, 
their thermal conductance is thus negligible and the electrical power dissipated in the resistor flows to the 
substrate through the unique gold wire (25 ^m in diameter, 3 mm long) ultrasonically bonded between 
the resistor and the gold pad evaporated on the sapphire substrate (on the right in Fig. 1). By measuring 
the temperature of the sample when a given power is injected in the resistor, this system gives an over- 
estimation of the thermal resistance between the sapphire substrate and the cryostat. It has the advantage 
of being reversible in comparison with a resistor deposited on the sapphire substrate, as the gold wire can 
be easily removed. 

We show now that heating phenomena are irrelevant to account for the I- V non linearities. We first 
focus on the possible heating of the whole sample together with its sapphire substrate, and then turn to 
"electronic heating" . 



The measurement of the thermal resistance TZth between the cryostat on one hand and the sample 
with the sapphire substrate on the other hand was carried out as follows. A constant voltage V = Fl 
was applied to the Y^Six-a sample, leading to a current I(F,T) depending on the cryostat temperature 
T. A given electrical power 1 nW <V< 100 nW was then dissipated in the small 1 Mf2 resistor. This 
power flowed to the heat sink through the substrate whose temperature was then increased by ST, leading 
to an increase SI of the current. We checked that SI was proportional to V. Assuming for a while that 
the TV non linearities do not result from sample heating, we extracted, from the data of Fig. 2, ST 
from 81 and deduced TZth = ST/V. Fig. 9 gives the resulting TZth a s a function of T. From it, we get 
7?.th(74mK) ~ lOmK/nW, while we can see on Fig. 2 that the I(F, 74mK) and I(F, 124 mK) curves begin 
to merge for IV ~ 0.04 nW which corresponds to ST ~ 0.4 mK much lower than the 50 mK separating 
these two 1(F) characteristics. This confirms the above assumption of irrelevance of sample heating; and 
it can be easily shown to be true at any temperature. Let us note that TZth exhibits a TZth ~ T~ 3 behavior 
which characterizes the Kapitza thermal resistance at the boundary between two materials. In our case, 
they are the sapphire substrate and the gold pad thermally connected to the cryostat. 



In this section, we show that the non-linearities of our I- V curves cannot be ascribed to a heating effect 
as found by Wang et al. for NTD sensors [Q. Assuming the validity of such a model for our sample, the 
power P = IV = IFl injected in the YSi "electron bath" would increase its temperature to a value T e 
larger than the cryostat temperature T, leading to a YSi electrical resistance 



A. Heating of the whole sample and substrate 



B. Electronic heating 



R(T, F) = i?o exp ((T /T e y) 



(Al) 



with T e given by 



P + P P =IV = g(T? - T n ) 



(A2) 
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where g and r\ are parameters characterizing the thermal conductance between the electron bath and the 
cryostat, and P v is the parasitic power injected in the sample due to the limited rf shielding etc. As we have 
excluded a possible heating effect of the sapphire substrate (see previous section) , the thermal resistance to 
consider is either the Kapitza one at the boundary between the YSi sample and the sapphire, or the thermal 
resistance due to electron-phonon coupling in the YSi itself. In the first case rj = 4p3^, while in the second 
one, 77 = 5 to 6 [0 J|^-[j53|] (64j. We investigated the experimental values of dP/dT e , which would be equal 
to grfT^~ l , hence would not depend on T and P p . Within the heating model, T e can be extracted for each 
(I,V) point using Eq. (Al). Fig. 10 shows dP/dT e as a function of T e . The fact that most of the curves 
depend on T is a strong argument against the heating model. However, we can see a trend towards a T 
independence for low T values. As indicated by the straight lines corresponding to power laws, these low T 
curves are not compatible with the heating model because they correspond to 77 values which depend on T e 
and may be much larger than 4 to 6. Finally, we note that in materials close to a-YSi, the electron-phonon 
coupling parameter g was found experimentally to be of the order of 10 3 W K -6 cm -3 (in NbSi, with rj = 6) 
[ p2[ or 200 W K~ 5 cm~ 3 (in AuGe, with 77 = 5) |34|. As a consequence, for our 7.2 x 10~ 5 cm 3 sample, 
we expect dP/dT e ~ 0.4 W/K or dP/dT e ~ 0.07 T e 4 W/K which is several orders of magnitude larger 
than what we find (see Fig. 10). If we assume that the possible heating effect is due to the YSi-sapphirc 
Kapitza resistance, we can use the experimental g ~ (1 — 10) x 10~ 3 W K~ 3 cm~ 2 values for sapphire-metal 
interfaces [^l), which lead to dP/dT e values of the order of 10 -3 T 3 W/K, again several orders of magnitude 
larger than our experimental values. Yet, those very large discrepancies guarantee that the heating effects 
are negligible. In comparison to data from authors who see heating effects, this can be explained by : i) 
our very resistive sample, ii) the low electron-phonon and Kapitza resistances due to the geometry of the 
sample. 
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APPENDIX B: ESTIMATION OF £ AND e r FROM T AND T. 



The detailed calculations of A, B, a, al and (3 were performed by Bottger et al. assuming Mott VRH 
(7 = l/(d + 1)) while we have 7 = 1/2. Thus our fit result A£ = 1.1 nm cannot be used straightforwardly 
to get £. In this appendix, we first try to estimate in which interval A must lie in the case of 7 = 1/2 and 
then we use the value of To as an additional experimental constraint on £. 

It appears that the exponents a, a' and /? should remain unchanged when going from 7 = l/(d+ 1) to 
7 = 1/2. This can be readily seen for a — a' — 1 in the framework of Pollak's calculations [||. Moreover, 
the fact that a — 1 for 7 = 1/2 is confirmed by many analyses of experimental results using Eq. (4) j2i| 
|jT| p5| fist [ JlfH pl| | J42] | (4f| ; while, in what concerns al the experimental situation is unclear (see Sec. 
I). However, the prefactors A and B are likely to be changed by the presence of a Coulomb gap. In this 
case, we can state at least A < 1/2, since A = 1/2 comes from the replacement of E m / (ksT) ~ ^(Tq/T) 7 
by E m / '(fc^T) — eFr m / '(fceT) in Eq. (1). Such a substitution obviously overestimates / in the presence 
of an electric field since it neglects both the influence of returns and the insensitivity of a large number of 
pair currents with respect to F. From A < 1/2 we deduce £ > 2.2 nm. Using the Bottger, Bryksin et al. 
value H @ A =1/6, we find £ = 6.6 nm which, using Eq. (5) gives e r = 29, while £ > 2.2 nm leads 
to e r < 86. 

According to Imry et al. |32|, one expects e r — Ni(£/ Xt.f.) 2 where Xt.f. is the Thomas Fermi screening 
length and Ni — 1 according to Abrahams and Lee |p5[ . Thus, using Eq. (5), At.f. — 0.08A -3 / 2 nm which 
leads to At.f. = 1.2nm for A = 1/6 and At.f. > 0.23nm for A < 1/2. In crystalline metals with a 
concentration n of one electron per atom we have typically At.f — 0.06 nm. Here we expect n to be of 
the order of 0.1 and since the standard screening theory yields At.f. oc n -1 / 6 , the value of 0.23 nm for our 
sample is plausible. However, we cannot exclude that the upper value of 1.2 nm is plausible too. Indeed, it 
was theoretically found |66]] that close to the transition, electronic diffusion is considerably lowered, which 
should reflect in a lowering of screening, i.e. in an (unknown) increase of Xt.f.- Moreover, if as recently 
suggested [|67|, many body effects come into play in Efros-Shklovskii VRH, one expect both a reduction of 
the 2.9 factor in Eq. (5) and a change of the A value. Finally, considering all these unknown effects, At; 
is found to have the correct order of magnitude and we estimate that £ ~ 7 nm. 

The determination of the dielectric constant e r of the system is a very interesting working direction 
for the future due to the fundamental interest of e r in localization and MIT studies. This is also a strong 
argument in favor of a precise theoretical determination of the numerical parameters in Eq. (3). Even 
if those parameters are not completely known, the relative evolutions of e r as a function of the dopant 
concentration, magnetic field, etc. should be reachable with the e r extraction method we used. 
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Figure Captions 



FIG. 1. Experimental setup: the 9 /im thick Yo.19Sio.8i layer is deposited onto a sapphire substrate 
with two interdigitcd gold electrodes. Thermalization of the sample was ensured by 43 gold wires (4 
represented) bonded on the left gold pad. Upper right corner : small 1 MCI chip heater allowing thermal 
resistance measurements at low temperature (see Appendix A). 

FIG. 2. Current as a function of electric field at various temperatures below 1.3 K : for clarity only 18 
out of the 36 different curves were represented. Inset: R(F — > 0) vs. 1/T (dots), for temperatures above 
4 K. The dashed line corresponds to an activated law with a characteristic temperature of 40 K. 

FIG. 3. Logarithm of the resistance as a function the electric field to temperature ratio, for 14 tem- 
peratures among the 28 measured ones. 

FIG. 4. Normalized value of T -7 for 3 values of 7 as a function of the logarithm of the resistance at 
zero electric field. The normalized T~ 7 is (T~ 7 - T^p in ) / (T^p ax - T^ 7 n ) where Tuin and Tuax are the 
two extreme values of the temperature T selected for the plot. Inset : residue (i.e. measured minus fitted 
value) of the linear fit of \n[R(F — > 0)] as a function of T~ 7 , for 7 = 1/4 (open circles) and 7 = 1/4 (closed 
circles) . 

FIG. 5. Effective exponent of the current vs electric field law, as a function of the electric field, for 6 
temperatures. 

FIG. 6. The length L extracted from the slopes of the ln(V/i") vs F/T lines at F -> using Eq. (4) 
(dots), and the fitted L — L (To/T) 1 ^ 2 law (continuous line). 

FIG. 7. The *(F,T) function (given by (T/T0) 1 / 2 ]n(I/I E .s.), with I E .s. = (V/Ro) exp[-(T /T) 1 /2]) 
as a function of F/T for all the data except the very high field points. \f increases more slowly than F/T 
(the dashed line corresponds to ^ oc F/T) and its evolution can be fitted (solid line) using the Bottger 
et al. prediction * = T(F/T) - A(F/Tf with = 1.15 ± 0.1, T = (1.25 ± 0.5) x 10~ 5 KmV -1 and 
AV/5/r = 0.70 ± 0.07. Inset : same * vs. F/T dependence with linear scales showing that the A(F/T)@ 
term plays a role at low F/T values; the dashed line corresponding to $ oc F/T. 

FIG. 8. F/T dependence of T(F/T) - (F, T) , with T = 1.25 x lO^KmV -1 , for all the data points 
except the very high field ones. The dashed line corresponds to a A(F/T) /3 dependence with (3 = 1, and 
shows that (i > 1. Solid line : fit with j3 = 1.15 and A = 1.4 x 10 _6 (Km/V)' 3 . Inset : same data with 
linear scales, the dashed line corresponding to a ~ F/T dependence. 

FIG. 9. Thermal resistance lZ t h(T) measured by dissipating a controlled power in the small 1 MO chip 
(see Fig. 1). Dashed line : fit of the data with a T~ 3 law. 

FIG. 10. Derivative of the power dissipated in the YSi sample with respect to the effective electron 
temperature T e , as a function of T e . For each (I, V) point, T e is calculated using Eq. (Al). The points line 
up in curves corresponding to different cryostat temperatures among which 5 are indicated in the figure. 
For comparison, the three straight lines correspond to 3 power laws dP/dT e oc Tj 7 " 1 ^ with 7 — 1 = 3, 5 
and 9. 
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